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Abstract 

It is demonstrated that local quark-hadron duality prescription applied to several exclusive 
processes involving the pion, is equivalent to using an effective qq (two-body) light-cone wave 
function SSf( LD '(x,k±) for the pion. This wave function models soft dynamics of all higher 
qG . . . Gq Fock components of the standard light-cone approach. Contributions corresponding 
to higher Fock components in a hard regime appear in this approach as radiative corrections 
and are suppressed by powers of a s /ir. 
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1 Introduction 



Study of effects due to finite sizes of the hadrons and incorporation of the transverse momentum 
degrees of freedom is a notoriously difficult problem for the QCD analysis of inclusive (see, e.g., 
0) and exclusive processes. Since the advent of the parton model 0, it is taken for granted 
that the hadron can be viewed as a collection of quarks and gluons, each of which carries a finite 
fraction XiP of the large "longitudinal" momentum P of the hadron, and also some "transverse" 
momentum ki±. However, a justification of such a picture from the basic principles of QCD is not a 
straightforward exercise. To begin with, the coordinate-representation version of k± is a derivative 
dj_ in the transverse direction. In a gauge theory, dj_ always comes together with the gauge field 
A± in the form of a covariant derivative d± — > D± = d± — igA±, i.e., the finite-size effects are 
mixed with those due to extra gluons. 

In the operator product expansion (OPE) approach ||, the nonperturbative aspects of the 
hadron dynamics are described/parameterized by matrix elements of local operators. In particular, 
the longitudinal momentum distribution is related to the lowest-twist composite operators in which 
all the covariant derivatives appear in traceless-symmetric combinations {D^D^ . . . -D Mn } @, ^ ||. 
To take into account transverse-momentum effects, one may wish to consider matrix elements of 
higher-twist composite operators in which some of the covariant derivatives appear in a contracted 
form like D 2 = D^JD^ . Attempting to relate them to transverse-momentum distributions (cf. [[7|), 
one would immediately notice, however, that D 2 looks more like analogue of the quark virtuality 
k 2 . Furthermore, the presence of the gluonic field in the covariant derivative obscures such an 
interpretation as well. In particular, using the equation of motion ^D^q = 0, one can convert 
a two-body quark-antiquark operator g{7 Ml D M2 . . . D fln }D 2 q with extra D 2 into the "three-body" 
operator <?{7 Ml -D M2 • • • D Mn }(cr Ml 'G ^ u )q with an extra gluonic field G^ u . Moreover, the two contracted 
covariant derivatives _D M . . . can be separated by _D Mfc 's forming the traceless combination and, 
to put _D M and _D M next to each other, one should perform commutations, each producing a G-field 
again 0. By choosing an optimized basis, one can avoid some of the complications |l|, ||, but the 
observations listed above clearly show that the OPE-inspired approach is unlikely to produce a 
simple and intuitively appealing basis for constructing phenomenological functions describing the 
transverse-momentum degrees of freedom. 

Still, the OPE approach has many evident bonuses. In particular, it is based on a covariant per- 
turbation theory in 4 dimensions and provides an explicitly gauge-invariant and Lorentz-covariant 
description. In this respect, it is analogous to the Bethe-Salpeter formalism for bound states. 
However, the well-known problem of the Bethe-Salpeter formalism is the presence of the unphys- 
ical variable of the relative time. This variable is unnecessary and, without a loss of information, 
one can describe bound states in a 3-dimensional formalism. This is achieved by projecting the 
Bethe-Salpeter amplitude on a particular (e.g., equal-time or light-front) hyperplane. 

In the light-cone approach (see, e.g., |10j), a hadron is described by a set of light-cone wave 
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functions (Fock components) ty( N '({xi, fcj^})- For mesons, the two-body wave function ^^(x, k±) 
can be related to the Bethe-Salpeter amplitude taken in the light-cone gauge and integrated over 
the minus component of the relative momentum |TTJ. A bonus of the light-cone approach is that 
effects due to transverse momenta are unambigously separated from those related to higher Fock 
components. Hence, if the lowest Fock component gives the dominant contribution, one can hope 
to construct a reasonable phenomenology based on modelling the two-body LC wave function. On 
the other hand, if one should perform a summation over all Fock components, the predictive power 
of the scheme is rather limited, since constructing models for numerous higher Fock components 
leaves too much freedom for the model builders. This poses a serious problem for phenomenological 
applications of the standard light-cone formalism. In particular, if one uses the popular BHL set 
of constraints for the pion LC wave function [12|, the qq component always contributes less than 
50% to the pion form factor value at Q 2 = 0, and higher Fock components are absolutely needed 
to ensure the correct normalization F 1T {Q 2 = 0) = 1. 

A possible way out is to introduce an effective two-body wave function (see, e.g., [IB]), which 
includes the low-energy contribution from the higher Fock components, so that one would get 
F-niQ 2 = 0) ~ 1 just from the overlap of these wave functions. One can interpret such a wave 
function as a wave function for a constituent quark, i.e., a quark dressed by soft gluons. However, 
for different processes, the higher Fock components can appear with different process-dependent 
weights, and it is not clear a priori whether the effective wave function can be introduced in a 
universal way. Another point is that while absorbing information about soft gluons into a qq wave 
function may be a good approximation, hard gluons cannot be correctly taken into account in this 
way. Hence, even using the effective two-body wave function, one should allow for a possibility 
of having explicit multi-body wave functions. Still, the dominant role of the effective two-body 
component may take place in such a scheme, since each emission of a hard gluon is suppressed by 
the QCD coupling constant a s /n ~ 0.1, and the contribution of the multi-body components may 
be relatively small. The problem is that it is unclear how to combine the QCD corrections with 
the constituent quark picture, because the constituent quark is not a field one readily finds in the 
original QCD Lagrangian. 

Here, using the pion as an example, we will outline a new approach to transverse-momentum 
effects in exclusive processes. It is based on QCD sum rule ideas. On several examples, we show 
that results obtained using the quark-hadron duality prescription |14] can be reformulated in terms 
of a universal effective wave function ty LD (x, k±) absorbing information about soft dynamics. The 
scheme starts with diagrams of ordinary covariant perturbation theory and allows for a systematic 
inclusion of the radiative corrections in a way totally consistent with the basics of QCD. 
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2 Handbag diagram and ^-scaling 



A naive idea is that, to take into account effects due to the finite size of the hadrons, one should just 
write the "parton model" formulas without neglecting intrinsic transverse momentum in hard scat- 
tering amplitudes. In doing this, however, one should explicitly specify a field-theoretic approach 
which is used for such a generalization of the standard parton model. As emphasized in the Intro- 
duction, one can choose here between at least two basically different alternatives: standard covariant 
4-dimensional formalism or 3-dimensional approaches analogous to the old-fashioned perturbation 
theory. The bonus of the 4-dimensional approach, in the form of the OPE, is a gauge-invariant 
and a Lorentz-covariant description of the hadrons in terms of matrix elements of composite oper- 
ators. However, as argued above, interpretation of the OPE results in terms of transverse degrees 
of freedom is not self-evident. Moreover, there are some practically important amplitudes which 
are "protected" from the dynamical (-D 2 )"-type higher-twist corrections. The most well-known 
example is given by the classic "handbag" diagram for deep inelastic scattering. As we will see 
below, in a scalar toy model its contribution contains only target-mass corrections, i.e., it gives no 
information about finite-size effects. In QCD, the handbag contribution contains a twist-4 operator 
with extra D 2 , but no operators with higher powers of D 2 . 

2.1 Scalar model 

To illustrate the effect in its cleanest form, let us consider the handbag contribution in a model 
where all fields are scalar (Fig.|l|): 

At large Q 2 = —q 2 , one can neglect the parton virtuality k 2 in {k + q) 2 = q 2 + 2{kq) + k 2 and 
expand the propagator in powers of 2{kq)/q 2 to obtain 

T(q,p)= fd'kf:^^F(p,k). (2.2) 

J n=0 ^ 

Now, the integral 

J ... k^F(p, k)d A k = A n p^ ...p^ + traces (2.3) 

is evidently the matrix element of a local operator with n derivatives. The usual parton density f(x) 
is introduced by treating the coefficients A n as its moments ||, |5| : 

A n = [ x n f(x)dx. (2.4) 
Jo 

As a result, the amplitude can be written as 
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a) b) c) 



Figure 1: Handbag diagram, a) Momentum representation for scalar model, b) Coordinate 
representation for scalar model, c) QCD modification of the quark propagator. 

f{x)dx + 0{l/Q' i ) 

Taking its imaginary part W(q,p) ~ImT(g,p), we get 

W(q,p) = f Q f(x) 5(Q 2 - 2x(qp)) dx + 0(1/Q 4 ) = ^x B f(x B ) + 0(1/Q 4 ), (2.6) 
where x B = Q 2 /2(qp) is the standard Bjorken variable. 

2.2 Power corrections and ^-scaling 

Eq.(|2.6p gives the lowest-twist contribution. The power-suppressed terms denoted by 0(1/ Q^) 
are apparently due to the neglected k 2 term in the original propagator. One can expect that 
supplementing the {kq)/Q 2 expansion by the k 2 /Q 2 expansion, one can take into account the 
effects due to nonzero virtuality k 2 by introducing phenomenological functions related to matrix 
elements like 

J k^ 1 . . . (k 2 ) N F(p, k)d 4 k = A^pf 11 traces . (2.7) 

Of course, one should be more careful now with the "traces" in this parameterization. The best way 
to maintain the necessary accuracy is well-known: one should take the traceless part . . . k^"} of 
the original tensor k^ 1 . . . k^ n . Then, the right-hand-side will also be a traceless tensor, constructed 
from the 4- vector p^ . 

So, if one decides to keep the k 2 terms, one should supplement this by a re-expansion of the 
(/cg) ri -factors over traceless tensors. In fact, the actual problem is simpler than it seems, because a 



T(q,p) 



1 f 1 00 



2x(qp) 
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straightforward expansion of the propagator is just in terms of the traceless combinations: 



(q + kY 



on 



k<q „=o (Q 



Note, that there are no {k^ . . . k fln }(k 2 ) terms with N 7^ in this expansion: the (k 2 ) N terms 
from a naive expansion over powers of (kq) and k 2 are exactly cancelled by {k 2 ) N terms from 
the reexpansion of (qk) "-factors over traceless tensors. In other words, the handbag diagram is 
insensitive to nonzero-virtuality effects. Introducing the twist- 2 distribution function via 

J{k^ ... }F(p, k)d 4 k = {p^ 1 ...p^yj x n f{x)dx (2.9) 

and performing the summation over n by inverting the expansion formula (|2.8f) , we get 



T(q,p)= / . - ,J (x)dx. (2.10) 
Jo [q + xpy 

The essential point is that no power-suppressed terms were neglected in this derivation. Hence, we 
can write (q + xp) 2 = —Q 2 + 2x(qp) + x 2 p 2 keeping all the terms here and calculate the imaginary 
part: 

W(q,p) = f f(x) 5 (-Q 2 + 2x(qp) + x 2 p 2 ) dx = X * /(0, (2.11) 

JO V If , 4p x B 



'1 + 

where 



2xr . 

(2.12) 



ip 2 x 2 



D_ 



is the Nachtmann-Georgi-Politzer £- variable [Jq], Hence, all the power-suppressed contributions 
contained in the handbag diagram can be interpreted as the target-mass corrections. In particular, 
the handbag contribution contains no power corrections for a massless target. 



2.3 Coordinate representaion 

Absence of the higher twist terms as well as the possibility to easily calculate the target-mass 
dependence of the handbag contribution is directly related to the fact that the propagator of a 
massless particle has a simple singularity structure. To illustrate this, let us write T(q,p) in the 
coordinate representation: 

T(q,p)= f^e^(p\m^)\p)- (2-13) 
The first term in the z 2 -expansion for the matrix element 

(pM0)V(*)b) = 6(zp) + z 2 U*P) + {z 2 YUzp) + ■■■ (2-14) 
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corresponds to the twist-2 distribution amplitude 



V (x)e ix{zp) dx 



(2.15) 



while subsequent terms correspond to operators containing an increasing number of <9 2 's. It is 
straightforward to observe that, while the twist-2 term produces the l/Q 2 contribution, the twist-4 
term is accompanied by an extra z 2 -factor which completely kills the 1 / ^-singularity of the quark 
propagator, and the result of the d A z integration is proportional in this case to 5 4 (q — xp), i.e., this 
term is invisible for large Q 2 . The same is evidently true for all the terms accompanied by higher 
powers of z 2 . This means that the handbag diagram contains only one term: it cannot generate 
higher powers of l/Q 2 which one could interpret as the {{k 2 )/Q 2 ) n expansion. 

For spin-1/2 particles, the quark propagator S c (z) ~ ^^/(z 2 ) 2 has a stronger singularity 
for z 2 = 0, which is cancelled only by the 0(z A ) term in the expansion of the matrix element 
(p\(j{0) , ~f t j,q(z)\p). Hence, one may expect that there is a non-vanishing twist-4 contribution corre- 
sponding to the 0(z 2 ) term of this expansion, but no higher terms. In a gauge theory, like QCD, 
one should also take into account the fact that the gluonic field A u , in a covariant gauge, has zero 
twist. As a result, if the gluons have longitudinal polarization, the configurations shown in Fig.|l]c 
are not power-suppressed compared to the original handbag contribution. The net result of such 
gluonic insertions into the quark propagator is a phase factor 

S c (z) -> S\z)Pe i9zV ti A » {tz)dt {l + 0(G)}, (2.16) 

where the 0(G) term corresponds to insertion of physical gluons. The latter are described by the 
gluonic field-strength tensor G a p and produce l/Q 2 -suppressed contributions. Thus, including the 
phase factor, we get the modified QCD handbag contribution 

T(q,p) ~ J d'ze-^ z ^{p^^^ ^ A ^ tz)dt q(z)\p). (2.17) 

The matrix element 

(p\q(0h,e l9zV Jo A '^ dt q(z)\p) (2.18) 

can be Taylor-expanded just like (p\q(0)'y^q(z)\p), with the only change d u — ► D v in the resulting 
local operators. Thus, the incorporation of gauge invariance does not change our conclusion that 
the (generalized) handbag diagram cannot generate a tower of the (1/Q 2 ) ra corrections which one 
could interpret as the ({k 2 ± )/Q 2 ) n or ((k 2 )/Q 2 ) n expansion. The power corrections are produced 
by the final-state interaction which is described by complicated contributions due to the operators 
of qG . . . Gg-type. At twist 4, the qD 2 q and qGq terms combined together can be interpreted 
in terms of functions related to operators of q~D\q type [|IJ. However, for twist-6 and higher, 
the absence of the generic q(D 2 ) N q contribution stops further progress in this direction. Hence, 
a simple phenomenological description of higher-twist corrections in terms of something like the 
transverse-momentum distribution f(x, kr) is impossible. 
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Figure 2: Form factor of the 7*7* — > 7r° transition, a) General structure, b) Leading-order pQCD 
term. 

3 Exclusive processes: 7*7* — ► tt° transition 

The transition 7*(<?i)7*(<72) — > n°(p) of two virtual photons into a neutral pion (Fig.||a) is the 
cleanest exclusive process for testing QCD predictions. The relevant form factor F^^^o (q 2 ,Q 2 ), 
with q 2 = —q 2 ,Q 2 = —q 2 , can be defined in terms of the pion-to-vacuum matrix element of the 
product of two electromagnetic currents: 

4tt J d 4 xe-^ x (n,P \T {J^x) J„(0)} |0) = iV2e m F rr7T o (q 2 ,Q 2 ) . (3.1) 



3.1 pQCD results 

In the lowest order of perturbative QCD, the asymptotic behaviour of R,» 7 * 7r o (q 2 , Q 2 ) can be 



calculated from a diagram similar to the handbag diagram (see Fig.|2|6) |T(J. The basic change is 
that one should use now the pion distribution amplitude ip n (x) instead of the parton density f(x): 



2 , Air f 1 ip w {x) 



3 Jo xQ 2 + xq 2 



rfx + 0(a s /7r) + 0(l/Q 4 ). 



(3.2) 



Experimentally, the most important situation is when one of the photons is real: q 2 = 0. In this 
case, pQCD predicts that [|Kj 



f^(Q 2 ) = y Si ^ dx + 0{asM + 0(1/g4) ' 



xQ 

The nonperturbative information is accumulated here by the integral 



(3.3) 



dx 



o x 



(3.4) 
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Its value depends on the shape of the pion distribution amplitude ip n (x). In particular, using the 
asymptotic form |L6], [17], [UJ 

rt s (x) = 6f n x(l-x) (3.5) 
gives the following prediction for the large-Q 2 behaviour [|T0[] : 

F^AQ 2 ) = (3-6) 
3.2 Anomaly and BL-interpolation 

Of course, the asymptotic l/Q 2 -dependence cannot be the true behaviour of F 77 * 7r o(Q 2 ) in the 
low-Q 2 region, since the Q 2 = limit of i 7 77 * 7r o(Q 2 ) is known to be finite and normalized by the 



7T° — ► 77 decay rate. In fact, incorporating PCAC and ABJ anomaly fl8 |, one can calculate 
F 77 *„-o(0) theoretically: 

FyyAty = —f~- ( 3 - 7 ) 
^ Jit 

In pQCD, one can imagine that the transition from the high-Q 2 asymptotics to the low-Q 2 behaviour 
is reflected by higher twist corrections of (M 2 /Q 2 ) n -type, which may sum up into something like 
1/(Q 2 + M 2 ), i.e., some expression finite at Q 2 = and behaving like 1/Q 2 for large Q 2 . This idea 



was originally formulated by Brodsky and Lepage jnj who proposed the interpolation formula 

F 77 * 7r o((5 2 ) = 7- Ttr^i ( 3 - 8 ) 



7T/ W (l + 



4^ 

which reproduces both the Q 2 = value (3.7) and the high-Q 2 behaviour ( |3.6| ) with the normal- 



ization corresponding to the asymptotic distribution amplitude ( |3.5|) . 
The BL-interpolation formula (|3.8|) has a monopole form 

F ir AQ 2 ) ~ V(i + Q 2 M)) 

with the scale so = 47r 2 / 2 w 0.67 GeV^ 2 , which is numerically close to the p-meson mass squared: 
m 2 ~ 0.6 GeV 2 . So, the BL-interpolation suggests a form similar to that based on the VMD 
expectation F 77 * 7r o(<5 2 ) ~ 1/(1 + Q 2 /m 2 ). In the VMD-approach, the p-meson mass m p serves as 
a parameter determining the pion charge radius, and one can expect that the tower of {sq/Q 2 ) n - 
corrections suggested by the BL-interpolation formula can be attributed to the intrinsic transverse 
momentum. 

3.3 Light-cone formalism and power corrections 

As noted earlier, the relative weight and interpretation of power corrections depends on a particular 
formalism used for a beyond-the-leading-twist extension of pQCD formulas. In the operator product 
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expansion approach, the lowest-order (in a s ) "handbag" contribution to the 77* — > 7r° form factor 
again has only pion-mass corrections: 

(Q 2 ) = ^ /* , 2 da. (3.9) 

77 v 7 3 Jo xQ^ + x(l — xjm^ 

This means that, in the OPE approach, the ("so"/Q 2 ) 7V -type corrections can come only from the 
qG . . . Gq operators, for which simple phenomenology is impossible. 

As an alternative to the covariant perturbation theory and OPE, one can use the light-cone 
(LC) formalism [10], in which effects due to the intrinsic transverse momentum kj_ are described 
by the light-cone wave function *&(x, k±). The LC formula for the 77* — > tt° form factor looks like 

(e± x q±)F^AQ 2 ) ~ / ~^ )] dxd 2 k^ , (3.10) 

where q± is a two-dimensional vector in the transverse plane satisfying q\ = Q 2 and e± is a vector 
orthogonal to q± and also lying in the transverse plane ||10|| . In the LC formalism, quark and gluon 
fields are on-shell. However, the invariant mass yVl of an intermediate state does not coincide with 
the pion mass. In particular, for the gg-component, M 2 = (k\ + m 2 )/x(l — x). Hence, integrating 
over k± is equivalent to integrating over the invariant masses (or LC-energies) of intermediate 
states, with ty(x,k±) specifying the probability amplitude for each particular "mass". 

Unfortunately, eq.( |3.10| ) also has little chances of producing the series of the (kj_) / Q 2 corrections, 
because the expansion 



1 1 00 



(xq± - k_i_) 2 x 



I Am I < xQ) + 6(\k ± \ > xQ) 



...qfik? ...hp } (3.1i; 



(xQ 2 ) n+1 {kl/x) n+1 

contains only traceless combinations {k± . . . k 1 ^ 1 }. Substituting it into eq.( 3.1C ) and using that the 
wave function *&(x, k±) depends on k± only through k\, we obtain that all terms of this expansion 
(proportional to Legendre's polynomials) vanish after the angular integration, except for the n = 
and n = 1 terms. As a result, the leading 1/Q 2 term is corrected only by the term resulting 
from integration of ^(x, k±)/kj_ over the region \k±\ > xQ: there is no tower of {(k 2 _)/Q 2 ) N 
power corrections in the handbag contribution. Thus, one is forced again to explain the (l/Q 2 ) N - 
corrections by contributions from higher qG . . . Gq Fock components, which unavoidably leads to a 
complicated phenomenology. 

The last but not the least comment is that the LC formalism is based on bound-state equations 
like — K ® ty, with K being an "interaction kernel". However, it is not clear whether such an 
equation has any justification in QCD outside perturbation theory. Moreover, it is known that 
QCD has a lot of nonperturbative effects: complicated vacuum, quark and gluon condensates, etc., 
which play a dominant role in determining the properties of QCD bound states. So, it is very 
desirable to develop a QCD description of hadrons in terms of functions similar to the bound state 
wave functions *$?(x, k±), but without assuming existence of bound state equations. Below, we 
outline our attempt to derive such a description from QCD sum rules and quark-hadron duality. 
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4 Basics of quark-hadron duality 



4.1 Outline of the QCD sum rule calculation of f n 

The basic idea of the QCD sum rule approach |20| is the quark-hadron duality, i.e., the possibility 
to describe one and the same object in terms of either quarks or hadrons. To calculate we 
consider the p^p^-part of the correlator of two axial currents: 

n""(p) = i J e^(0\T (j+ (z)& (0) )| 0) d'x = p,p„U 2 (p 2 ) - g^ip 2 ). (4.1) 
The dispersion relation 

I r°° o( s) 

U 2 (p 2 ) = - / -t-^-t-rds + " subtractions" (4.2) 

IV JO S — p 2 

represents H^Qo 2 ) as an integral over hadronic spectrum with the spectral density p hadron (s) deter- 
mined by projections 

(0|^(0)|7r;P) = i/ ff P M , (4.3) 
etc., of the axial current onto hadronic states 

p hadron {s) = 7cf 2 5(s - m 2 n ) + Ttf 2 Al 8{s -m 2 Ai ) + "higher states" (4.4) 

(f° xp ~ 130.7 MeV in our normalization). On the other hand, when the probing virtuality p 2 is 
negative and large, one can use the operator product expansion 

U 2 (p 2 ) = IirV) + ~M GG ) + ^s(qq) 2 + ■■■ (4-5) 

where H 2 Krt {p 2 ) = li^^ip 2 ) is the perturbative version of Il 2 (p 2 ) given by a sum of pQCD Feynman 
diagrams while the condensate terms (GG), (qq), etc., (with perturbatively calculable coefficients 
A, B,etc. ) describe/parameterize the nontrivial structure of the QCD vacuum. 

For the quark amplitude Y[ < ^ Lark (p 2 ) , one can also write down the dispersion relation ( |4.2j ), with 
p(s) substituted by its perturbative analogue p quark (s): 

p quark (s) = -^(l + — + ..) (4.6) 

47T \ 7T / 

(we neglect quark masses). Hence, for large — p 2 , one can write 

1 roc p hadron^_ quark ^ A B 2 

-/ 2 —ds=-(a s GG) + -a s (qq} 2 + .... 4.7 

TT Jo S — p z p 4 p° 

This expression essentially states that the condensate terms describe the difference between the 
quark and hadron spectra. At this point, using the known values of the condensates, one can try 
to construct a model for the hadronic spectrum. 
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In the axial-current channel, one has an infinitely narrow pion peak p n = 7rf 2 5(s — m 2 ), a rather 
wide peak at s ~ 1.7 GeV 2 corresponding to A\ and then "continuum" at higher energies. The 
simplest model is to treat A\ also as a part of the continuum, i.e., to use the model 

phadron^ ^ _ ^ + p auark ^ ^ > ^ (4 _g } 

in which all the higher resonances including the A\ are approximated by the quark spectral density 
starting at some effective threshold Sq. Neglecting the pion mass and requiring the best agreement 
between the two sides of the resulting sum rule 

h = -\ - — 2 ^ s + -M GG ) + -Mw) 2 + ■■■ (4.9) 

p 2 n Jo s — p 2 p A p b 

in the region of large p 2 , we can fit the remaining parameters /„- and so which specify the model 
spectrum. In practice, the more convenient SVZ-borelized version of this sum rule (multiplied by 

M 2 ) 

is used for actual fitting. Using the standard values for the condensates (GG), (qq) 2 , we adjust s 
to get an (almost) constant result for the rhs of eq. ([4.1Cl|) starting with the minimal possible value 
of the SVZ-Borel parameter M 2 . The magnitude of f n extracted in this way, is very close to its 
experimental value f^ xp ~ 130 MeV. 



4.2 Local duality 

Of course, changing the values of the condensates, one would get the best stability for a different 
magnitude of the effective threshold so, and the resulting value of /„- would also change. There 
exist an evident correlation between the values of /„- and Sq since, in the M 2 — > oo limit, the sum 
rule reduces to the local duality relation 

fl = - H p* uark ( s )ds. (4.11) 



7T JO 

Thus, the local quark-hadron duality relation states that, despite their absolutely different 
form, the two densities pi uark ^s) and p /lctdron ( s ) give the same result if one integrates them over 
the appropriate duality interval sq. The role of the condensates was to determine the size of the 
duality interval Sq, but after it was fixed, one can write down the relation ( [4.1 1|) which does not 
involve the condensates. 

Using the explicit lowest-order expression pl uark (s) = 1/47T, we get 

s = 47c 2 f 2 . (4.12) 
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Notice that Sq = 47r 2 / 2 is exactly the combination which appeared in the Brodsky-Lepage interpo- 
lation formula ( |3.8| ). Numerically, 47T 2 / 2 « 0.67 GeV 2 , i.e., the pion duality interval is very close 



to the p-meson mass: m 2 « 0.6 GeV 2 . In fact, in the next-to-leading order 

P^) = ^(l + v)' ( 4 - 13 ) 

So, using a s /ir « 0.1, one gets so practically coinciding with m 2 . For the form factors, this leads 
to results close to the VMD expectations, even though no explicit reference to the existence of the 
p-meson is made. 

4.3 Local duality and pion wave function 

In the lowest order, the perturbative spectral density is given by the Cutkosky-cut quark loop 
integral (see Fig.|3]a) 



p quark {s) = — 



/^(l-^) (**) ((? " k ?) d ' k ( 4 - 14 ) 



2tt 2 

where s = p 2 . Introducing the light-cone variables for p and k: 

p= {p+ = P,p^ = s/P,p± = 0} ; k = {k + = xP, 
and integrating over we get 

pi uark (s) = f 1 dx [ 5 ( s - % J d 2 £; ± . (4.15) 



27T 2 Jo 7 V XX 

The delta-function here expresses the fact that, since we are working in the 4-dimensional formalism, 
the light-cone combination k\jxx coincides with s, the square of the external momentum p. 
Substituting now p? uarfc ( s ) into the local duality formula, we obtain 

fl = ^£ dx J 6 ( k l^ x * s °) Sk ^ ■ ( 4 - 16 ) 
This representation has the structure similar to the expression for f n in the light-cone formalism 



U = VQ Cdx [*( x ,k ± )^±, (4.17) 

JO J o7T 

where \&(x, k±) is the gg-component of the pion light-cone wave function. To cast the local duality 
result (|4.16j ) into the form of eq. (|4.17| ), we introduce the "local duality" wave function for the pion: 



V LD (x, k ± ) = ^ 0{k\ < xxs ) . (4.18) 

Jit 
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a) 



b) 



Figure 3: Leading-order contributions for a) two-point spectral density p(s) (4.14) and b) three- 
point spectral density p(s\, $2,Q 2 ) (5.8). The narrow dashed lines indicate Cutkosky cuts. 

The specific form dictated by the local duality implies that ty LD (x, k±) simply imposes a sharp 
cut-off at k\xx = so. 

It should be emphasized that in eq.( [4.17D we are integrating over k±, i.e., the combination k\xx 



no longer coincides with the mass 2 of the external particle (which is m 2 in our case). This is 
precisely the feature of 3-dimensional bound-state formalisms: internal particles are on mass shell 
k 2 = m 2 , but off the energy shell: k\xx ^ m 2 . In what follows, we show, on less trivial examples, 
that the local duality prescription allows one to get results reminiscent of 3-dimensional formalisms, 
though all the actual calculations are performed in the standard 4-dimensional perturbation theory. 



5 Pion electromagnetic form factor 
5.1 Sum rule 

To demonstrate that the function ty LD (x, k±) really has the properties one expects from the pion 
wave function, let us consider the quark-hadron duality in a more complicated context of the pion 
electromagnetic form factor F n (Q 2 ). It is defined by 

(p 2 |J r/i (0)bi) = (rf+rf)F ir (g a ), (5.1) 

where Q 2 = — Q»2 ~Pi) 2 - To apply the QCD sum rule technique, we should consider in this case 
the correlator [14]] 

Zj(pi,ftO = i J e^^{0\T{j p {y)J^(0)j+{x)}\0)d i xd i y (5.2) 

of two axial currents and one electromagnetic current J M . The pion EM form factor can be 

extracted from the invariant amplitude T{p\,p\ ) Q 2 ) corresponding to the structure P^P^P^, where 
P = (pi+p 2 )/2. 

14 



The obvious complication now is that we have two channels to be "hadronized" , since the pion 
is present both in the initial and final states. This necessitates the use of the double dispersion 
relation 

T(»?,»o.fl 2 ) = — / dsi [ ds 2 7 — ^ Sl ow ' — - + " subtractions" (5.3) 

7l 2 Jo JO (Si - p{){s 2 ~ Pi) 

involving the double spectral density p(s\, s 2 ,Q 2 ). Its hadronic version p hadron (si, s 2 , Q 2 ) contains 
the term corresponding to the pion form factor 

82, Q 2 ) = 7T 2 f 2 F 7T (Q 2 )5(s 1 - ml)5(s 2 - m*) (5.4) 

and the contributions corresponding to transitions between the pion and higher resonances, and 
also the terms related to elastic and transition form factors of the higher resonances. To construct 
the two-dimensional analog of the "lowest state plus continuum" ansatz, we will treat all the 
contributions, except for the p T7T , as "continuum", i.e., we will model p hadron (si, s 2 ,Q 2 ) by the 
p quark (si, s 2 ,Q 2 ) outside the square (0, s ) x (0,s ): 

p(si, 8 2 , Q 2 ) = p™{si, s 2 , Q 2 ) + (1 - 9(8! < s )6(s 2 < s )) /^'(si, s 2 , Q 2 ). (5.5) 
The SVZ-borelized sum rule (with M\ = Mf = M 2 ) for the pion form factor then has the form 



21, 14 



1 r s o r&o 

f 2 F n (Q 2 ) = — d Sl ds 2 p^ k { Sl ,s 2 ,Q 2 )^ V 

7T Z JO JO 



S_! + S 2 

M 2 



+ 12^ + 81^^ 13+ M^ (5 ' 6) 



(the pion mass was neglected as usual). 



5.2 Local duality 

In the large-M 2 limit, this gives the local duality relation [[14 



fX D {Q 2 ) = - 2 Hds! r ds 2 p^ k {s!,s 2 ,Q 2 ). (5.7) 

7T Z JO JO 

Again, the perturbative spectral density p quark (s!, s 2 , Q 2 ) corresponding to the triangle diagram 
Fig.^ can be easily calculated using the Cutkosky rules and light-cone variables in the frame where 
the initial momentum pi has no transverse components pi = {p^ = P,Pi = si/P, 0j_}, while the 
momentum transfer q = p 2 — pi has no "plus" component: p 2 = {P, (s 2 + Q 2 ±)/P, Q±}'- 



^\s x , s 2 , Q 2 ) = A jf 1 dx J d 2 k ± 5 L - 5 (s 2 - {k±+ J q? ) ■ 



(5.8) 
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Here, x is the fraction of the total "plus" light-cone momentum carried by the quark absorbing the 
momentum transfer from the virtual photon and k± is its transverse momentum. 

Substituting p quark (s 1 , s 2 , Q 2 ) into the local duality relation, we get the light-cone formula for 
the pion form factor 

F^ D (Q 2 ) = £ dx * LD {x, fcj V LD (x, k ± + xq), (5.9) 



where ty LD (x, k±) is exactly the local duality ( [4.18| ) wave function introduced in the previous 
section. 

At Q 2 = 0, the LC formula reduces to the integral 

which can be interpreted as the probability to find the pion in the state described by the wave 
function ty(x, kj_). Using the explicit form of ^> LD (x, k±), we immediately get 



J 7T 



which reduces to 1 if we take the lowest-order value s = 4n 2 f 2 for the duality interval. Hence, 
p(°) = 1, i.e., the probability to find the pion in the state described by *ff LD (x,k±) is 100%. In 
other words, ^ LD (x, k±) may be treated as an effective wave function absorbing the low-energy 
information about all Fock components. 

It should be emphasized, however, that if one uses the next-to-leading order value 



s d) _ ^fi 

~~ 1 _|_ K 

for the duality interval, the probability integral P will be smaller than 1 (?W ~ 0.9 for a s ~ 0.3). 
This is a direct manifestation that the local duality prescription explicitly produces contributions 
which can be interpreted as hard parts of the higher Fock components like qGq, etc. (see Fig.|| and 
Fig.[5] below). However, the total probability to find the pion in such a higher Fock state is rather 
small: it is suppressed by the factor — sa 0.1 . 

Though the probability integral P differs from 1 beyond the leading order, the relation ^^(0) = 
1 holds to all orders of perturbation theory. The reason is that, at any order, there exists the 
Ward identity relation between the 3-point function T^g(p,p) and the 2-point function Ii a p{j)): 
T£g(p,p) = — dTlapip) / dp^. As a result, the 3-point spectral density p guark (si, s 2 , Q 2 ) reduces to 
the 2-point spectral density p quark (s): 



P 



quark (s u 8 2 , Q 2 = 0) = n8( Sl - s 2 )p quark ( Sl ). (5.11) 
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Hence, for Q 2 = 0, the duality integral for the pion form factor automatically reduces to a one- 
dimensional integral over s\ (or s 2 ), and the duality integral (|5.7| ) for f 2 F n (0) would coincide with 
that for f 2 ( |4.11| ), provided that the sides of the duality square for the three-point function are 
exactly equal to the duality interval Sq for the two-point function (the latter not necessarily being 
equal to the lowest-order value Sq° = An 2 / 2 ). As a result, the local duality prescription gives 
F£ D (0) = 1 to all orders of perturbation theory. In the lowest order, it also gives = 1. 

One should not overestimate the accuracy of the local duality results in the region of small Q 2 . 
Though F^ D (Q 2 ) dictates the values rather close in magnitude to the VMD curve F^ MD {Q 2 ) or 
any other fit to data, the LD-formula 



F^ D {Q 2 ) = 1 - 




(5.12) 



gives infinite slope at Q 2 = 0, and one should not use it for calculating the derivatives of F n (Q 2 ) 
below Q 2 ~ So- As emphasized above, we obtained the correct value for F n (0) only because this 
value was protected by the Ward identity. It is well known that the 3-point function in the small-Q 2 
region has more complicated quark-hadron duality properties which require a separate study. For 
the same reason, the local duality fails to produce reasonable valence parton densities for the pion. 



6 Quark-hadron duality for the F] ' * (Q 2 ) form factor 
6.1 Basics 

Within the QCD sum rule approach, one can extract information about the 7*7* — > 7r° form factor 
from the three-point correlation function (||22||): 

Air r 

^(gi,g 2 ) = d 4 xd 4 ye-^ x -^(0\T{J^x)J u (y)j 5a (0)}\0) (6.1) 

calculated in the region where all the virtualities q\ = ~q 2 ,q\ = —Q 2 and p 2 = (qi + q 2 ) 2 are 
spacelike. 

To study the form factor F 7 * 7 *„.o(g 2 , Q 2 ), one should consider the invariant amplitude F (p 2 , q 2 , Q 2 ) 
corresponding to the tensor structure efj, uqiq2 p a . The dispersion relation for the three-point ampli- 
tude 

F (V, g 2 , Q 2 ) = - r P(g,g2,< p 2) ds + subtractions (6.2) 

V ' 71 JO S — p 2 

specifies the relevant spectral density p(s,q 2 ,Q 2 ). Its hadronic version 

p hadron ^ ^ q 2^ = ^^.^(^ q^ S (s - m\) + "higher states" (6.3) 
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contains the term with the form factor we are interested in. The relevant perturbative spectral 
density p quark (s,q 2 ,Q 2 ), in the lowest order, is given by the integral representation 



p^ rk (s, q\ Q 2 ) = 2 f 1 5(3 - q2xiX3 + Q2x ' X A 5(l-j2x i )dx 1 dx 2 dx 3 
Jo \ x\Xi ) i=1 



(6.4) 



in terms of the Feynman parameters for the one-loop triangle diagram. Scaling the integration 
variables: x\ + x 2 = y, x 2 = xy, x\ — (1 — x)y = xy and taking trivial integrals over x 3 and y, we 
get 

p^\s, g 2 , Q 2 ) = 2 f 1 x ^f 2 +^ ") 2 dx . (6 .5) 
Jo [sxx + xQ 1 + xq 2 \ 6 

It can be shown that the variable x here is the light-cone fraction of the pion momentum p carried 
by one of the quarks. 

6.2 Local duality 

Incorporating the local duality, we obtain 

i^%(g 2 , Q 2 ) = ^ f P QUark ^ « 2 > Q 2 ) = J r[ ldx r ds r ^ ( ! Q ntf 2)2 213 ■ ( 6 - 6 ) 
' ' Txj-K Jo hJk Jo Jo [sxx + xQ z + xq z \ A 

Substituting the variable s (the mass 2 of the qq pair) by the light-cone combination k^/xx, we 
get F^P^ (q 2 , Q 2 ) as an integral over the longitudinal momentum fraction x and the transverse 
momentum k±: 



Finally, introducing the effective wave function ^ LD (x, kj_) given by (|4.18|) we can write F LD (Q 2 ) 
in the "light-cone" form: 

It is instructive to analyze this expression in some particular limits. 
6.3 Limits 

1. Both photons are real. When both q 2 and Q 2 are small, we can use the fact that 



// 1 



(/1 2 + k\f 2 



-5(k{) (6.9) 
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in the /i 2 — > limit, to obtain (cf.[[L^]) that the 7r° — > 77 decay rate is determined by the magnitude 
of the pion wave function at zero transverse momentum: 

F 7 ^o(0,0) = — L= £ ^ LD {xM = 0)dx. (6.10) 

Using the explicit form of \l/ LD (x, k±) ( |4.18| ), we obtain 

F 7 ^(0,0) = -L (6.11) 

which is exactly the value ( |3.7| ) dictated by the axial anomaly. 

2. pQCD limit. Assuming that both q 2 and Q 2 are so large that the /c^-term can be neglected, 
we get the expression 

F - M , Q 2 ) = -J— J* _i£_ J f kx ^ t j + 0(1/Q <) . (6,2) 

Identifying the wave function integrated over the transverse momentum with the pion distribution 
amplitude 

ri D (x) = -0- J * LD (x, k ± ) d 2 k ± , (6.13) 
we arrive at the lowest-order pQCD formula 

for the large-virtuality behaviour of the 7*7* — > vr° transition form factor. 

3. One real photon. A very simple result for p quark (s, q 2 , Q 2 ) appears when q 2 = 0: 

p^ fc ( S ,g 2 = 0,Q 2 ) = ^^. (6.15) 

This formula explicitly shows that if Q 2 also tends to zero, the spectral density p quark (s,q 2 = 
0,Q 2 ) becomes narrower and higher, approaching S(s) in the Q 2 — > limit (cf. [p3fl ). Thus, the 
perturbative triangle diagram dictates that two real photons can produce only a single massless 
pseudoscalar state: there are no other states in the spectrum of final hadrons. As Q 2 increases, the 
spectral function broadens, i.e., higher states can also be produced. Assuming the local duality, we 
obtain: s 

hF^AQ 2 ) = - f P quark ^ °> ^ ds = mxw v (6 ' 16) 

" vr Jo 7r(l + Q z /s J 

For large Q 2 , this gives 

F^AQ 2 ) = ^ + 0(l/Q 4 ). (6.17) 
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This result can also be obtained from the q 2 = version 

4tt f 1 <p v (x) 



f p ^{Q 2 ) = yJ ^r dx + °(VQ 4 ), ( 6 - 18 ) 



xQ< 

of the pQCD formula ( |6.14j ), if we will use there the asymptotic form of the pion distribution 
amplitude 

^ D (x)=Qf 7T x(l-x) (6.19) 

produced by the local duality prescription. 

In other words, the local duality formula (|6.16 ) exactly reproduces the Brodsky-Lepage inter- 



polation ( |3.8| ) between the Q 2 = value 1 / 7rf n fixed by the AB J anomaly and the leading large-Q 2 
term 4irf n /Q 2 calculated for the asymptotic form of the pion distribution amplitude. 

7 Normalization properties of the effective wave function 

7.1 Momentum representation 

By construction, ty hn ( ) satisfies the standard constraint 

imposed on the two-body Fock component of the pion light-cone wave function by the correspon- 
dence with the 7r — ► fiu rate. 

Furthermore, the x-integral of ty LD (x, k±) at zero transverse momentum 

f 1 dx^ LD (x,k ± = 0) = (7.2) 
Jo U 

has the right magnitude to produce the correct value F 7 . 7 * 7r o(0, 0) = l/nf^ imposed by the 7r° — > 77 
rate. Note, that this value is by a factor of 2 larger than the constraint imposed in |T2|| on the quark- 



antiquark component of the LC pion wave function. The difference can be traced to the claim, 
made in jl2fl , that the gg-component of their pion wave function gives only a half of the 7r° — > 77 
decay amplitude. The other half, it was argued, should be attributed to the gg7-component of the 
pion wave function. Within our approach, an analogue of the gg7-component appears only in the 
first order in a s . At the leading order, there is only one term, which, I repeat, correctly reproduces 
the Q 2 = value of F 7 . 7 «„.o(0, Q 2 ). 



7.2 Probability integral 

The integral 



P = I^^M*,^)\ 2 . (7.3) 
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gives the probability to find the pion in the state described by the wave function *&(x, k±). As dis- 
cussed earlier, substituting the local duality wave function ^ LD (x, k±) into this relation, one would 
get P = s /(47r 2 f 2 ) = 1. At the lowest order in a s , the local duality wave function ty LD (x, k±), in 
this sense, describes 100% of the pion content, i.e, ^ LD (x,k±) can be understood as an effective 
wave function dual to all Fock components of the pion light-cone wave function. 

It is worth noting here that, in our approach, P < 1 for any wave function k±) which 

a) depends only on the combination k\jxx: ^(x, kj_) = f(k\/xx), 

b) monotonically decreases with the increase of k\/xx, 

c) never becomes negative and 

d) satisfies the constraints ( |7.1|) and (\l.2\). 

The upper limit for P is reached when k±) assumes the steplike form (|4.18| ) dictated by the 



local duality. The requirement that the (generalized) valence content of the pion should not exceed 
100% is not unreasonable. Furthermore, I fail to see why, in a particular model, this probability 
cannot reach 100%. However, if instead of our constraint ( |7.2| ) one apllies that proposed in [JT2]] , the 
upper limit for P, under the same conditions a) — d), is 0.5, i.e., one should mandatorily require 
that at least 50% of the pion content must always be related to non-valence components. 

7.3 Impact parameter representation 

Defining the impact parameter b± as the variable which is Fourier-conjugate to the transverse 
momentum k±: 

■$(x,k ± ) = ( e ik ^^(x,b ± )d 2 b ± , (7.4) 



we can write down the normalization conditions for the 6j_-space wave function ty(x,b±). Eq. 
(|7. 1|) , following from the requirement that n — ► \iv rate is specified by gives the magnitude of 
\I/ LD (x, b±) at the origin: 

[ l dx^ LD (x,b ± = 0) = ^, (7.5) 
Jo V6 

and eq. ( |7.2j ) , following from the requirement that ir° — > 77 rate is given by axial anomaly, specifies 
its integral over the whole 6^-plane: 

1 dx f V LD {x, & l) d 2 b ± = (7.6) 

J u 

There is a widespread opinion that the axial anomaly is a purely short- distance phenomenon 
produced by ultraviolet divergences. However, the constraint (|7.6|) involving integration over all 
impact parameters clearly shows that the axial anomaly is deeply related to the long-distance 
physics as well. In particular, calculating the spectral density p(s,q 2 ,Q 2 ) exhibiting the anomaly 
behaviour in the q 2 , Q 2 — > limit, we never faced any ultraviolet divergences (cf. |23fl). 
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a) b) c) d) 



Figure 4: Some two-loop contributions to the two-point spectral function p(s) (8.1). a,b) Con- 
tributions corresponding to 0(a s ) corrections to the two-body qq effective wave function. c,d) 
Contributions corresponding to presence of hard three-body gGg-components both in the initial 
and final states. 

For reference purposes, we also give the impact-parameter version of the formula for the pion 
electromagnetic form factor: 

F^ D (Q 2 ) = -L f 1 dx ( e ix ( Q±b ^\$ LD (x,b ± )\ 2 d 2 b ± (7.7) 

4:71 JO J 

and the 6j_-space form of our effective wave function: 

^ LD (x,b ± ) = ——-y/xXSoJ^y/xXSo), (7.8) 

where J\{z) is the Bessel function. 



8 Higher-order corrections 

Calculating the spectral densities p quark (s, . . .) to higher orders in a s , we can study effects due to 
gluon radiation. Depending on the position of Cutkosky cuts, one can interpret, e.g., the next- 
to-leading order contributions either as corrections to the two-body qq effective wave function 
(Fig.|]a, b) or as three-body qGq Fock components (Figfjc, d) . 

In practice, even the lowest 0(a s ) correction requires a two-loop calculation, which is rather 
involved, especially for three-point functions. For the two-point function, the correction is known 
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Perots) = P Q Uark (s) + P r rk (s) = i- (l + ^) . (8.1) 

According to the Ward identity ( |5.11| ), this result can be used to get the Q 2 = value of the 
0(a s ) contribution to the spectral density p quark (si, S2, Q 2 ) related to the pion electromagnetic 
form factor. As a result, the 0(a s )-correction to the pion form factor for Q 2 = is given by 

5F K) ( g2 = 0) = ^(fo). ( 8 .2) 

7T 
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a) b) c) 



Figure 5: Some two-loop contributions for the spectral function p(si, S2,Q 2 )- a) Correction to 
the electromagnetic vertex, b) Contribution corresponding to the three-body gGg-component in 
the initial state and two-body gg-component in the final state, c) Term producing the 0(1/ Q 2 ) 
contribution. 



The duality interval so, in this case, is a natural (and the only possible) scale for the running 
coupling constant. 

Another important piece of information can be obtained for large Q 2 . In this limit, in con- 
trast to the one-loop term Po Uark (sx, s 2 , Q 2 ), which decreases like l/Q 4 , the two-loop contribution 
Pi Uark (si, s 2 , Q 2 ) contains a term (Fig.[5|c) which behaves like l/Q 2 - 



P T aTk (s u s 2 , Q 2 ) = ^p q uark (s 1 ) P r rk (s 2 ) + 0{l/Q% (8.3) 

where p < Q Lark {s\) and pl uark (s 2 ) are the lowest-order two-point function spectral densities (see 
eq. ( |3 . 1| ) ) . This behaviour agrees with the pQCD factorization theorem fT6| , |T0]] and quark counting 
rules. Substituting the asymptotic expression for pf iark {si, s 2 , Q 2 ) into the local duality relation 
(|5.7|) , we get the large-Q 2 behaviour of 5F^. a ^(Q 2 ): 



i.4) 



This result corresponds to the pQCD formula for the one-gluon-exchange contribution to the pion 
form factor R pi, HU 



Ff CD (Q 2 ) = [dxf^dy ^(x) vAv)^- 2 , (8-5) 

if one uses the "asymptotic" distribution amplitude <p® s (x) = 6f n x(l — x) dictated by the local 
duality. Now, by analogy with the Brodsky-Lepage interpolation, we can construct a model for 
5F^ as \Q 2 ) based on the simplest interpolation between its Q 2 = value and large-Q 2 asymptotics: 

5F^(Q 2 ) = (*) - ■ -L- , (8.6) 



7T 



1 + Q 2 /2s 
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Combining the 0(1) and 0(a s ) terms, we get the next-to-leading order LD-model for the pion form 

™ = (8.7) 

7T 

where F^ D ^ (Q 2 ) is the lowest-order result given by eq.( 5.12|) . For a s /7r one can take a constant 
value a s /n = a s (s )/n ~ 0.1 though, for truly asymptotic Q 2 , the scale of a s should have a 
Q 2 -dependent component. The curve based on eq. (|8.7|) is in good agreement with existing data. 



9 Conclusions 

Our main goal here was to demonstrate that the results of the approach based on local quark- 
hadron duality and QCD sum rules can be reformulated in terms of the effective wave function 
ty LD (x, k±) describing both longitudinal and transverse momentum distribution of quarks inside 
the pion. This approach has the following features: 

1) It is directly related to the QCD Lagrangian, and all the calculations are based on Feynman 
diagrams of standard 4-dimensional perturbation theory. 

2) As a result, the approach is fully compatible with high-Q 2 pQCD calculations and other QCD 
constraints, e.g., those imposed by the axial anomaly. 

3) Radiative (higher-order in a s ) corrections can be added in a regular way, through a well-defined 
procedure. 

4) There is no need for a special procedure separating soft vs. hard terms. In a sense, they are 
separated automatically by the duality interval parameter s - The "hard" terms have a natural 
subasymptotic modification in the low-Q 2 region. 

5) The bulk (soft) part of the higher-twist effects is described by an effective 2-body wave function 
fy LD (x, k±) rather than by increasingly complex wave functions of higher Fock components. 

6) In this approach, the contributions which can be interpreted as effective wave functions for the 
higher Fock components are small because they are suppressed by powers of a s (so)/vr. Hence, the 
effective valence component dominates and, in this respect, this approach resembles the constituent 
quark model. However, there is no need to introduce constituent quark masses. The scale respon- 
sible for the IR cut-offs is set by the duality interval s . 

7) The effective wave functions are introduced in this approach without any appeal to the existence 
of bound state equations. 

Local quark-hadron duality was also applied to nucleon form factors [031 and to 7^ — ► A 
transition form factors [J25| . The results of these studies can be used to develop a similar formalism 
for the baryons. Another possible development is to substitute the steplike effective wave functions 
\I/ LD (x, k±) by smooth functions, but without violating the constraints ( |7. 1|) , (|7.2|) and = 1. 
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